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Abstract 

We present further analysis of an anisotropic, non-singular early universe model that leads to the viable 
cosmology presented in JTJ. Although this model (the DLH model) contains scalar field matter, it is rem- 
iniscent of the Taub-NUT vacuum solution in that it has biaxial Bianchi IX geometry and its evolution 
exhibits a dimensionality reduction at a quasi-regular singularity that one can identify with the big-bang. 
We show that the DLH and Taub-NUT metrics are related by a coordinate transformation, in which the 
DLH time coordinate plays the role of conformal time for Taub-NUT Since both models continue through 
the big-bang, the coordinate transformation can become multivalued. In particular, in mapping from DLH 
to Taub-NUT, the Taub-NUT time can take only positive values. We present explicit maps between the 
DLH and Taub-NUT models, with and without a scalar field. In the vacuum DLH model, we find a periodic 
solution expressible in terms of elliptic integrals; this periodicity is broken in a natural manner as a scalar 
field is gradually introduced to recover the original DLH model. Mapping the vacuum solution over to 
Taub-NUT coordinates, recovers the standard (non-periodic) Taub-NUT solution in the Taub region, where 
Taub-NUT time takes positive values, but does not exhibit the two NUT regions known in the standard 
Taub-NUT solution. Conversely, mapping the complete Taub-NUT solution to the DLH case reveals that 
the NUT regions correspond to imaginary time and space in DLH coordinates. We show that many of the 
well-known 'pathologies' of the Taub-NUT solution arise because the traditional coordinates are connected 
by a multivalued transformation to the physically more meaningful DLH coordinates. In particular, the 
'open-to-closed-to-open' transition and the Taub and NUT regions of the (Lorentzian) Taub-NUT model 
are replaced by a closed pancaking universe with spacelike homogeneous sections at all times. 
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I. INTRODUCTION 

In a previous work [0Q we argued that it is natural to consider a generalisation of the standard 
cosmological scenario to one in which a scalar field dominates the dynamics of a homogeneous 
but, in general, anisotropic (Bianchi) universe. We presented a new solution (the DLH model) to 
the cosmological field equations based on a closed biaxial Bianchi IX universe containing scalar 
field matter. This led to a nonsingular 'pancaking' model in which the spatial hypersurface vol- 
ume goes to zero instantaneously at the 'big-bang', but all physical quantities, such as curvature 
invariants and the matter energy density remain finite, and continue smoothly through the big- 
bang. Moreover, we showed that the model leads to a viable cosmology at late times, exhibiting 
desirable features such as isotropisation and inflation, as well as producing perturbation spectra 
consistent with observations. 

We also noted in [1J that, despite containing scalar field matter, our model was reminiscent of 
the Taub-NUT vacuum solution, since both have biaxial Bianchi IX geometry and an evolution that 
exhibits a dimensionality reduction at a quasi-regular singularity. In this paper, we show that the 
metrics for the DLH and Taub-NUT models are, in fact, related by a coordinate transformation. 
It is thus of interest to investigate the explicit mapping between models based on each metric, 
with and without scalar field matter. Moreover, we investigate the well-known 'pathologies' of 
the Taub-NUT solution, in the context of the mapping to the DLH model. We contend that the 
natural coordinatisation of the DLH model is more physical than the traditional coordinates used to 
describe Taub-NUT. We thus consider the possibility that pathologies such as the 'open-to-closed- 
to-open' transition and the Taub and NUT regions of the (Lorentzian) Taub-NUT model arise 
because the traditional coordinates are connected by a singular transformation to the physically 
more meaningful DLH coordinates. 

This paper is organised as follows. We begin with a description of Bianchi universes in Section 
|n} We then briefly review the DLH model in Section [HI] and the Taub-NUT model in Section IV 



We show that the two metrics are related by a coordinate transformation in Section |V} where the 
mapping naturally leads to a new reparameterised form of the Taub-NUT metric. We investigate 



this reparameterised Taub-NUT model in Section VI and find that it admits a simple scaling family 



of solutions, in contrast to the conventional Taub-NUT setup. In Section VII[ we consider the 



mapping between the DLH and reparameterised Taub-NUT models, with and without a scalar 



field, and interpret the solutions physically, before concluding in Section VIII 



II. BIANCHI MODELS 

Bianchi universes are spatially homogeneous and therefore have a 3 -dimensional group of 
isometries G3 acting simply transitively on spacelike hypersurfaces. The standard classification 
hence follows Bianchi's classification of 3-parameter Lie groups 0. 

We adopt the metric convention (H ). Roman letters a,b,c... from the beginning of the 

alphabet denote Lie algebra indices. Greek letters /i,v,<7... label spacetime indices, whereas 
Roman letters i,j,k... from the middle of the alphabet label purely spatial ones. 

The isometry group of a manifold is a Lie group G and can be thought of as infinite simally 
generated by the Killing vectors t,, which obey [<^u,£v] = C^ V ^ CT where the CZ V are the structure 
constants of G. These can be used to construct an invariant basis, which is often useful to make 
the symmetry manifest. This is a set of (basis) vector fields X^, each of which is invariant under 
G, i.e. has vanishing Lie derivative with respect to all the Killing vectors such that 

[S fl ,X v ]=0. (1) 

Such a basis can be constructed simply by imposing this relation at a point for some chosen set of 
independent vector fields and using the Killing vectors to drag them out across the manifold. The 
integrability condition for this set of first-order differential equations in fact amounts to demanding 
that the C?j v be the structure constants of some group. The invariant vector fields satisfy 

|Xjx,X v J = —C^ v X(j. (2) 

Denoting the duals of the X^ (the so-called invariant 1 -forms, or Maurer-Cartan forms) by 00^, the 
corresponding curl relations for the dual basis are 

dco^ = ]-C^ t (0 a A (0 T . (3) 

Because the X^ are invariant vectors, the metric can now be expressed as 

ds 2 = g^ lv CO^(O v 1 (4) 

for some g^y 

Bianchi models can be constructed in various different ways. For simplicity, we use the fact 
that the timelike vector generating the foliation of spacetime into homogeneous spacelike hyper- 
surfaces commutes with the three Killing vectors within the hypersurfaces (generating the homo- 
geneity), and hence we choose a representation that is diagonal: 

ds 2 = dt 2 - Yij(t)(O l (o j = dt 2 - y k i(t)(e k i (x)dx l )(e l j(x)dx j ) 1 (5) 



in terms of an explicit (x,y,z) coordinate system. 

As outlined in (H, the Bianchi classification of G3 group types hinges on the decomposition 
into irreducible parts of the spatial part of the structure constants C\-. Imposing the Jacobi identity 
then essentially leaves nine distinct choices of parameterisations (zeroes and signs) of the structure 
constants corresponding to nine different groups, called Bianchi I through Bianchi IX. All Bianchi 
models have a timelike vector, which generates the preferred foliation into spacelike hypersurfaces, 
and three spacelike Killing vectors, generating the homogeneity on those hypersurfaces. Both 
models that we will consider, DLH and Taub-NUT, are Bianchi IX such that the symmetry algebra 
is5o(3) ~su(2). 

General Bianchi IX models are thought, generically, to exhibit complicated dynamics and 
chaos. The dynamics near the initial singularity of vacuum and orthogonal perfect fluid mod- 
els is believed to be governed by Bianchi I and II vacuum states via the Kasner map ('Mixmaster 
attractor'). This description in terms of successive Kasner periods can be reformulated in terms of 
reflections and is called 'cosmological billiard motion', also known as 'BKL analysis' after Belin- 
skii, Khalatnikov and Lifshitz [3]-[9]|. This turns out to be just an example of a more general phe- 
nomenon when one considers (super-)gravity close to a spacelike singularity (the 'BKL- limit'). In 
this limit the gravitational theory can be recast in terms of billiard motion in a portion of hyperbolic 
space, as above, such that the dynamics is determined by successive reflections. These reflections, 
however, are precisely the elements of a Lorentzian Coxeter group, which are themselves the Weyl 
groups of (infinite-dimensional) Kac-Moody algebras. This then leads to the conjecture that these 
Kac-Moody algebras are in fact symmetries of the underlying gravitational theory |fT0] - PT2ll . There 
are also some concerns about the discrete nature of the Kasner map, and a continuous generalisa- 
tion - see, for instance [fT3]4T5l . For recent work on locally rotationally symmetric (LRS) Bianchi 
cosmologies with anisotropic matter see, for example, ifToll . In HI, we remarked on how Bianchi 
models can be considered as a deformation of Friedmann-Robertson- Walker (FRW) models, and 
observed how these perturbations freeze out during inflation. In fact, this can be understood in 
terms of the characterisation of Bianchi IX models as an FRW model deformed by long range 
gravitational waves [fT7|. 

For the DLH and Taub-NUT models that we will consider, however, there exists an additional 
biaxial symmetry: two of the left-invariant SU(2) one-forms appear with the same coefficient 
in the metric. Thus there is an additional right action by a U(l) factor inside the SU(2) which 
acts by isometries, so we are considering a class of metrics admitting an SU(2) x U(l) symmetry 



group. As shown in fl}, and demonstrated further below, this additional symmetry allows for much 
simpler dynamical evolution than in the full triaxial Bianchi IX case. 

III. THE DLH MODEL 

For the DLH model, the metric of the form (Q is 

ds 2 = dt 2 - \R\(co 1 ) 2 - \R\ [(co 2 ) 2 + (co 3 ) 2 ] , (6) 

which trivially reduces to FRW form in the special case where the two scale factors are equal, 
R{(t) = R 2 (t). Following flU, and with the usual definitions for the Hubble parameters Hj(t) = 
Ri/Ri for the different directions, the Einstein field equations are easily computed and give two 
dynamical equations for the two independent radii R\ and R2 (from now on we will drop the 
explicit t dependence of variables): 

2H 2 + 3H 2 + K P -A= 1(3^-4) (7) 

and 

2H 2 + 2H l ~H 2 + 2H l H 2 + Kp-A=-^(5^-4\ (8) 

as well as the Friedmann equation (or Hamiltonian constraint) 

H} + 2H l H 2 -K P -A=±(^-4), (9) 

and equation of motion for a simple massive scalar field (for which the potential is given by 

V(<j>) = \m 2 § 2 ) 

m 2 (j) + (H l +2H 2 )<i) + ^=0. (10) 

It turns out that these equations have relatively straightforward series solutions in t — to. We 
choose to = for simplicity, and take it to denote a big bang-like event. In 0]| we have presented 
two solutions with definite parity - one even (bouncing) and one odd (pancaking) in the non- 
degenerate scale factor R\. Here we are interested in the pancaking solution 

R x (t)=t (a + a 2 t 2 + a 4 t 4 + . . . ) 
R 2 (t)=R 3 (t)=bo + b2t 2 + b 4 t 4 + ... , (11) 

Ht)=fo + fit 2 + f4t 4 + ... 



ln(/? 2 (0) 



ln(/?,(0) 




FIG. 1: Dynamics of the DLH biaxial Bianchi IX model: evolution of the logarithm of the scale factors R\ 
and R2 in Planck lengths l p versus log time (t in units of Planck time t p ). 



where the dynamical equations ([7]), ([8]) and (10) allow one to fix the higher-order coefficients in the 
series order-by-order in terms of the initial values ciq = R\ (0) , bo — R2 (0) and /o = (0) . The fact 
that this also satisfies the Friedmann energy constraint ([9]) then proves that this odd-parity series 
solution is a valid expansion around the big-bang at t = 0, which one can use as a starting point 
for numerical integration. 

Fig. [T]shows the evolution of the scale factors R\ and R2 (for t > 0) for the viable cosmological 
solution presented in [Q, which is defined by the initial parameters ciq = 1.2, bo = 18000 and fo = 
13 (set by imposing a boundary condition at temporal infinity on the total elapsed conformal time, 
as suggested by lfl8l ). together with K = 1 and m = 1 /64000 (set in order to fix the normalisation of 
the resulting perturbation spectrum) and A set to zero, as it is dynamically unimportant at the early 
time scales that we are interested in. This cosmological model was obtained for a representative 
set of parameter values of order unity, rather than having been fixed in order to get best agreement 
with current data. These natural values were then scaled to the seemingly less natural values given 
above, using a scaling property discussed in the next paragraph. In order to fix the normalisation 
of the perturbation spectrum, the mass of the scalar field has to be rescaled and bo changes to a 
less natural value accordingly. This choice for the mass of the scalar field needs to be put in by 
hand for every current inflation model so does not constitute any unusual fine-tuning. 



As noted in (T), given a solution to the equations (J7])-(jT0]), a family of solutions is generated by 
scaling with a constant a and defining 

1 - - 9 

Ri{t) = —Ri(at), Hi{t) = aHi(at), ty(t) = <j>(at), m = am, A = or A. (12) 

LA/ 

This scaling property is valuable for numerical work, as a range of situations can be covered by 
a single numerical integration. Furthermore, many physically interesting quantities turn out to be 
invariant under changes in scale. This scaling property does not, however, survive quantisation, as 
any quantisation prescription for the scalar field introduces a length scale, which breaks the scale 
invariance. Therefore one would have to be careful when considering vacuum fluctuations. 

In [0Q| we demonstrated that at the time of pancaking, there is an instantaneous reduction in 
the number of dimensions of the homogeneous hypersurfaces, without a geometric singularity, or 
any singularities in physical quantities. Geodesies can extend through this point, though some of 
them may wind infinitely around the topologically closed dimension. In the light of the generic 
BKL-analysis mentioned above, it is interesting to note that the dynamics of our model is very 
straightforward. 

IV. THE TAUB-NUT MODEL 

The Taub-NUT model lfT9l is a biaxial Bianchi IX vacuum solution. Traditionally, the form of 
the metric chosen to describe it does not have the form ([5]), but is instead written 

ds 2 = 2dua) i -g(u)(a) l ) 2 -e 2 ^^[(a} 2 ) 2 + (Q} 3 ) 2 ], (13) 

where, for later convenience, we denote the standard Ryan & Shepley time coordinate by u (rather 
than t) and replace their function /3(f) by £(w). 

The metric ( fT3| ) can be shown to solve the vacuum Einstein equations provided 

. . Au + l-4B 2 u 2 C{u ) ( n 2 1 \* 

'( K > = *(W + i) ' e ={ Bu + Vb) ' (14) 

where A and B are arbitrary constants, which when varied lead to a family of solutions. The above 
functions are plotted in Fig. [2Jfor the choice of values A = B = 1. Note that e^ u > is always positive 
as asserted in lfT9l . but also that it behaves like \u\ for large values of B or u, whereas g(u) has 
the form of an inverted parabola for small u and approaches a constant for large u. Essentially, 
B measures the smoothness of e*( u > at the origin, whereas A shifts the centre of the approximate 




FIG. 2: Analytic Taub-NUT solution ( 14 1 for A = B = 1. 



parabola g(u). We note that there is no analogue of the simple scaling family of solutions ( 12) in 
this setup. 

The usual interpretation of Fig. [2] is that Taub-NUT has two NUT regions, corresponding to 
negative values of g and one Taub region, where g is positive. One infers that the Taub-NUT 
solution can be represented as a disc that evolves into an ellipsoid and back into a disc. In partic- 
ular, it is considered to evolve from timelike open sections in a NUT region, via lightlike sections 
(called Misner bridges), to spacelike closed sections in the Taub region, back into timelike open 
sections in the other NUT region. This open-to-closed-to-open transition is not mathematically 
singular, but it is incomplete, as geodesies spiral infinitely many times around the topologically 
closed spatial dimension as they approach the boundary [|20l - l26l (see, for instance, [27J for a recent 
treatment). This type of singularity is called 'quasiregular' in the Ellis and Schmidt classification 
||28l |29ll (these include the well-known 'conical' singularities [30]), as opposed to a (scalar or 
non-scalar) curvature singularity. The Taub-NUT solution therefore shows the same feature of di- 
mensional reduction as the DLH model and similarly does not have a geometric singularity during 
this collapse. Note, however, that the homogeneous hypersurfaces are only spacelike in the Taub 
region, leading to problems with Taub-NUT as a description of a homogeneous universe as one 
approaches the boundaries of the Taub region (the Misner bridges). 
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V. RELATIONSHIP BETWEEN DLH AND TAUB-NUT MODELS 

The similarities between the DLH and Taub-NUT models are sufficiently striking that the con- 
nection between the two models warrants further consideration. In particular, it is known that 
Taub-NUT is the only Petrov type D homogeneous closed vacuum spacetime and that all Petrov D 
solutions are known [31]. It is a straightforward, if tedious, calculation to show that the DLH met- 
ric ([6]) is also of Petrov type D. Moreover, one finds that the DLH and Taub-NUT metrics have the 
same degeneracy structure in their principal curvatures, i.e. the eigenvalues of the Riemann tensor 
(rather than the Weyl tensor used in the Petrov classification), in that there are two degenerate 
pairs and two singlets amongst the six real eigenvalues, as further explained in appendix |A} Fur- 
thermore, the most general biaxial Bianchi IX metric is also known to be the Plebahski-Demiahski 
metric ll3~2l . The similar behaviour of geodesies in both models is also suggestive. An analysis 
of geodesies in both models is contained in appendix [Bj and more details can be found in fl]. 
These similarities suggest that the DLH and Taub-NUT metrics might describe the same space- 
time geometry. In fact, as we now demonstrate, the two metrics are indeed related by a coordinate 
transformation. 

In general, the diffeomorphism invariance of general relativity means it can be difficult to de- 
termine if two metrics describe genuinely different spacetime geometries or are the same up to a 
diffeomorphism. Nonetheless, in the latter case, such a diffeomorphism can be found from con- 
sidering several scalar invariants. In general, four independent scalar invariants allow us to fix 
the diffeomorphism between the two spacetimes. A further invariant can then be used to check 
consistency or to derive a contradiction. In our present case, both models are homogeneous and 
so scalar invariants are functions of time only. Thus, by considering how a single curvature invari- 
ant, e.g. the Ricci scalar, transforms, one can straightforwardly identify an appropriate coordinate 
transformation linking the two metrics, which amounts simply to a time rescaling. 

Starting with the DLH metric (|6]), let us consider the time rescaling transformation 

u= f tl -Ri{t')dt' = r\t) (15) 

JO 2 

and define the functions P/(w) = Rj(f(u)) = Ri(t). The lower limit on the integral is chosen for 
convenience such that u = when t = 0. One could make another choice, but this simply adds a 
constant shift to the value of u. Moreover, let us shift the non-degenerate spatial one-form by a 
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timelike part to obtain 



a l = co l -\ --dt, o 2 = co 2 , o 3 = co 3 . (16) 

Ri{t) 



The shift is necessary to absorb unwanted terms arising from the time reparameterisation into a 
redefinition of the one-forms. As long as the one-forms are purely spatial to start with, the timelike 
component that we have added to the first one-form commutes through, such that this redefined 
set of one-forms still obeys the su(2) commutation relations and is therefore a valid set to describe 
homogeneous hypersurfaces. One thus obtains the metric 

ds 2 = 2duo l - \Pi (u) 2 (o 1 ) 2 - \P 2 (u) 2 [(o 2 ) 2 + (a 3 ) 2 ] . (17) 



Comparing this expression with the Taub-NUT metric (13), one sees that they have the same 



form, but with ^Pf replacing g and ^P| replacing e 2< = . It is straightforward to verify that the 



time rescaling (15) maps all curvature scalars for the DLH metric, such as the Ricci scalar (by 



construction), the Euler- Gauss -Bonnet invariant 

R GB= R iiva-c-^ R iiv J r R > ( 18 ) 

the Chern-Pontryagin scalar 

K C P = * R »vo^ va \ (19) 

and the eigenvalues of the Riemann and Weyl tensors (of course they are not necessarily indepen- 
dent) into those obtained for the Taub-NUT metric (c.f. appendix [A]). Moreover, we note that the 
form of the resulting Einstein field equations do not depend on the concrete realisation of the su(2) 
1-forms used in the metric; one simply requires that they obey the su(2) algebra commutation re- 
lations, which we have guaranteed by construction. A mapping between the geodesic equations in 
both models, for the explicit realisation used in 0], is exhibited in appendix |B} 



The time rescaling transformation (15) is valid independently of any concrete choice for the 
metric functions. Nonetheless, it can be seen that whenever R\[t) (or P{(u)) goes through zero, 
as is the case at pancaking events of the DLH model (or on the Misner bridges in the Taub- 
NUT model), this coordinate transformation will be problematic. The transformation itself is not 
singular, but one sees that when the integrand i?i changes sign (such as at the pancaking events), 
the definition of u becomes multivalued: the value of the integral for u will begin to decrease as t 
increases. Thus there is only a one-to-one correspondence between u and t as long as R\ and Pi 
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do not change sign. Furthermore, the definition of the 1 -forms (16) goes singular at the pancaking 
events also indicating a problem with u as a measure of time (c.f. Section VII D ). 
We also note that from ([15]) we have that 

du Ri(t) P\{u) 



dt 2 2 

and thus the inverse transformation is simply given by 

2 



(20) 



du'; (21) 



/o Pi(u') 

Hence DLH time t is essentially conformal time for Taub-NUT. This seems rather odd, as DLH is 
the natural generalisation of closed FRW models, and one is usually interested in conformal time 
associated with those cosmologically interesting solutions, which would make the usual cosmo- 
logical conformal time doubly conformal Taub-NUT time. 

VI. THE REPARAMETERISED TAUB-NUT MODEL 



The new form ( 17 ) for the Taub-NUT metric differs significantly from the traditional form ( 13 ). 
We believe that our new form in terms of the scale factors is physically more meaningful, as the 
squares of the scale factors premultiply the invariant one-forms, as in the DLH metric, or in the 
FRW special case, to give physical distances on the homogeneous spacelike slices. 



The differences between our form ( 17 ) and the traditional form ( 13 ) brings into question the 



usual interpretation of the Taub-NUT vacuum solution ( 14). In particular, we see that the tradi- 
tional g function is rather unnatural, as it corresponds to the square of a scale factor, rather than 
the more physically meaningful scale factor itself. Moreover, as the square of a real number, the 
function g should always be positive. Hence we should not allow the g function to go negative, 
and must instead select the positive branch at all times. This yields a very different picture from 
the alleged open-to-closed-to-open transition. If one simply took the modulus of the solution in 
Fig. [2j such that the ^-function is positive throughout, the resulting model is not a solution of 
the Einstein equations, unless e 2 ^ is also allowed to flip sign and become negative, thus raising 
a new problem. Rather, one should solve afresh the Einstein equations using the reparameterised 
Taub-NUT metric as the Ansatz. 



In our parametrisation ( 17 ), in terms of the scale factors Pi(u), with associated Hubble functions 



Kj(u) = Pi/Pi (where a prime denotes d/du) and scalar field F(u) = (/(«)) = (j>(t), the Einstein 

13 



field equations yield the dynamical equations 
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4P?K[ + SP 2 K 2 - 2P\k\ - 4km 2 F 2 + kP 2 F' 2 + 4P 2 K { K 2 -8A +—5^-4=0 (22) 

and 

8 / P 2 \ 
4P 2 K 2 + 6P 2 Kl - 4km 2 F 2 + kP 2 F' 2 + 4P\K y K 2 -8A -^3^-4=0, (23) 

°2 V °2 / 

as well as the Friedmann equation (or Hamiltonian constraint) 



yli/l i a-,. — „2p2 i ^d2z7/2 ^ d2 f. f. i oa i _ I _1 

A 2 V p 



2P{K{ + 4Km l F L + kP(F' z - 4PfK l K 2 + 8A +-^1-^-4)= 0, (24) 



2 



and, in general, the equation of motion for the scalar field 

P 2 F" + 2P 2 {K X + K 2 )F' + 4m 2 F = 0. (25) 

It is straightforward to show that, as expected, these equations also result from directly applying the 
time rescaling transformation ( p~5j ) and associated function redefinitions to the evolution equations 
©-([TO]) of the DLH model. 



We will solve the above system of equations, with and without scalar matter, in Section VII For 
the moment, however, we concentrate on the issue of scaling solutions. Considering the family of 
solutions in the DLH model related by { fT2] ), one could use the time rescaling transformation ( fT5| ) to 
map these solutions into our reparameterised Taub-NUT model, thereby constructing an analogous 
family of Taub-NUT solutions. In general, however, these will not be related by a simple scaling 



relation as in ( 12). Nonetheless, our reparameterised version of Taub-NUT does admit directly a 



family of solutions that are related by a straightforward scaling of the form 

Pi(u) = -— fl(j8»), K i (u)=pK i (pu), F(u)=F(Pu), m=^/pm, A = j8A. (26) 

VP 



Comparing this result with the corresponding scaling invariance (12) of solutions in the DLH 
model, we see that they both arise from a simple constant time rescaling of the form t —$■ at (in 
the DLH model) or u — > j5u (in the reparameterised Taub-NUT model). The relative square root 
between a and /3 results from the fact that the Taub-NUT metric is linear in the 'time' parameter 
u whereas the DLH metric is quadratic in t. 



The integral definition (15) of u(t) precludes a straightforward identification of u(ttt) with 



0Cu(t). One can, however, extend the diffeomorphism ( 15 ) to take us from scaled DLH to scaled 

Taub-NUT by writing 

rat i 

[3u= -Ri{t')dt' (27) 

Jo 2 
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and 



2a 

co 1 H ^—zdt, a 2 = co 2 , a 3 = co 3 , (28) 

Ri(at) 



and defining Pj( f3u) = Rt(at) andFQSw) = §(at). 

VII. COMPARISON OF DLH AND REPARAMATERISED TAUB-NUT MODELS 

Although there exists a coordinate transformation linking the DLH and reparameterised Taub- 
NUT metrics, this transformation is multivalued when pancaking events occur in the DLH model 
and at the Misner bridges in the Taub-NUT model; this leads to differences in the physical inter- 
pretation of the corresponding cosmological solutions. In this section, we therefore consider, in 
turn, the four cases of vacuum and scalar field matter solutions in both DLH and reparameterised 
Taub-NUT. We contend that the DLH set-up is the more physically meaningful. Anticipating this 
conclusion, we start by considering the DLH vacuum model as the most fundamental setup. 

A. DLH vacuum solution 



Setting the scalar field = in the DLH evolution equations (|7])-( 10 ), and also assuming A = 



for simplicity, as it is unimportant dynamically at early times, yields the system 

2H\ + 2H X - U\ + 2HiH 2 = ~ (5^ - 4j , (30) 

^ 2+2 ^ 2 =4Gf~ 4 )' (3i) 

which can be solved analytically, albeit in terms of elliptic integrals. The detailed derivation of the 
analytic form is given in appendix |Cj In short, the Einstein equations can be integrated to find one 
scale factor R{ and time t in terms of the other scale factor R2 as 

R 2{I )U» (32) 

and 



' Xl+a K+^E(yU) + \[^(ao + 2i)F(yU)- x [^(al-A)n(y;—^ ) , (33) 



2ao <2q + 4 V Si ' V 32 ' y 32ao \ ' iao 
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FIG. 3: Periodic 'DLH vacuum' solution for the values ciq = 1 and bo = 1. 



where E(z;k), F(z;k) and IT(z; V;&) are Legendre's three normal forms, /i an integration constant 
that can be fixed in terms of the initial conditions, and x and y are functions of i?2 defined in ap- 
pendix[C] Alternatively, the system of evolution equations can be solved numerically; comparison 



of the numerical and analytical solutions shows very good agreement (see Fig. 15 ). 

The solution is periodic in t, as shown in Fig. [3j which is a numerical solution with boundary 



conditions analogous to the pancaking solution in the case with the scalar field (11), i.e. 



*i(0) = 0, tfi(0) =ao, Ri(0) =b , # 2 (0) =0, 



(34) 



and we make the convenient choice ao = 1 , bo = 1 . In addition to the periodicity, the solution is 
symmetric about the two pancaking points in each period, with a parity inversion in R\ (which is 
linear near pancaking events) and Rj being even, in agreement with our previous pancaking DLH 
solution in (TJ. This is an interesting 'cyclic' model that repeats indefinitely. As we will see in 
the following section, it is stable to inclusion of a perturbing scalar field, but as the scalar field 
becomes heavier and/or denser, strict periodicity is broken. 
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B. DLH solution with a scalar field 

This physical set-up is, of course, that which we originally considered in 0}. As re-iterated 
in Section |ni} for some (quite natural) assumed values of the initial conditions and scalar field 
mass, such a model leads to a viable cosmology. It is of interest here, however, to investigate 
the transition from the cyclic DLH vacuum solution outlined above to the viable cosmological 
model by 'gradually' introducing the scalar field. This can be achieved by allowing the scalar field 
to become progressively denser or the mass of the scalar field heavier. Here the same boundary 



conditions are assumed as previously for the pancaking series solution (11 ), i.e. oddness for R\ 
and evenness for R 2 and <j>. 

In order to assess the effect of increasing scalar field energy density, the remaining parameters 
are kept constant at K = 1, A = 0, m = 1/64000, ciq = l,bo= 1, whilst 0o is varied over the range 
< 0o < 2 x 10 5 . Fig. [4] shows the vacuum solution in panel (a), and a small perturbation thereof, 
with 0o = 4, in panel (b). As the scalar field is increased, R[ turns around (panel c) and inflation is 
produced (panel d) (both for 0o = 1.22 x 10 5 ). Note that, from panel (d) onwards, we increase the 
range in t and take logarithms of the scale factors, to account for the fact that inflation is produced. 
Panels (e) and (f) show how higher initial scalar field energy densities produce more inflation (for 
00 = 1-3 x 10 5 and 0o = 2.0 x 10 5 respectively). 

Similarly, in order to study the effects of increasing the mass of the scalar field (Fig. [5]) (rather 
than its density), the other parameters are kept fixed at K = 1, A = 0, 0o = 10, ao = 1, ^0 = 1» 
whilst m takes the values: (a) m = 1, (b) m = 3, (c) m = 5, (d) m = 10 and (e) m = 100. Panel (f) 
shows the m = 100 case for a wider range in t. 

For a very light or diffuse scalar field, the equation of motion 

m 2 (j) + (H x + 2H 2 ) + 0=0 (35) 

is approximately satisfied by a constant scalar field, as the mass term is suppressed. This essen- 
tially eliminates the scalar field from the problem, so that we recover the vacuum case with its 
periodic behaviour. As the scalar field gets denser (see Fig. |4]) or the mass of the scalar field heav- 
ier (see Fig. [5]) (but still subject to the same boundary conditions), the deviations from the vacuum 
case grow, and eventually can no longer be considered small. The behaviour is then no longer 
periodic, and smoothly changes qualitatively, eventually yielding the inflationary cosmological 
solution described in 0}. 
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FIG. 4: Gradually introducing a light scalar field by increasing 0q, whilst keeping the other parameters fixed 
at K = 1, A = 0, m = 1/64000, oq = 1, bo = l: a) shows the vacuum model 0o = 0; b) 0o = 4; c) and d) 
both have <j) = 1.22 x 10 5 ; e) <j) Q = 1.3 x 10 5 and (f) 0o = 2 x 10 5 . 

C. Taub-NUT vacuum solution 



The usual derivation of the Taub-NUT vacuum solution is obtained by substituting the standard 



form of the metric ( 13 ) into the Einstein equations to yield the family of solutions ( 14 ). We now 
have two further methods for constructing this solution: by mapping the DLH vacuum solution 



using the diffeomorphism (15); and by substituting our reparameterised form (17) of the Taub 



NUT metric into the Einstein equations and solving the resulting system of equations given in 



Section VI These two alternatives approaches offer different physical insights into the nature of 
the Taub-NUT vacuum solution, and are considered below. 
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FIG. 5: Gradually increasing the mass of the scalar field whilst keeping the other parameters fixed at K = 1, 
A = 0, 0o = 1, oq = l,bo = l: (a) shows the model for m = 1; (b) m = 3; (c) m = 5; (d) m = 10; (e) and (f) 
are both for m = 100, with the latter covering a wider range in t. 

1. Mapping the DLH vacuum solution using the diffeomorphism 

Conceptually, the simplest way to arrive at the vacuum solution for the reparameterised Taub- 
NUT metric is to apply the diffeomorphism ( p~5| ) to the DLH vacuum solution discussed in Sec- 



tion VII A In practice, however, this rather complicated, as it involves elliptic integrals. For the 
sake of simplicity, we therefore concentrate on the period near (one of) the pancaking events, 
which will be sufficient to unearth some interesting properties of the mapping. 
In the vicinity of a pancaking event, R\ has the simple linear behaviour 



Ri(t)=a t. 



(36) 



Hence, the corresponding coefficient in the DLH metric (pL which is proportional to R\ °c t 2 , is 
a smooth function that touches zero at the pancaking (see Fig. [6]). To map this solution to the 



reparameterised Taub-NUT model, we first use the definition of u in the diffeomorphism ( 15 ) to 
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FIG. 6: Near-pancake limit of the DLH solution: R\{t) and R\(t) 
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FIG. 7: Near-pancake limit of the Taub-NUT solution: Pi (u) and Pf(u) 



obtain 



1 



1 



ciQt dt = -aot 
o 2 4 



(37) 



in this limit. Thus u is always positive (assuming ciq > 0), whilst t takes both positive and negative 
values. Also, since Ri(t) = P\(u), this yields 

pf(u)=alt 2 = 4aou. (38) 

We see immediately that the innocuous pancaking of the DLH solution now appears as ±y/u 
behaviour in the corresponding Taub-NUT solution (see Fig. [7]); this also matches the lowest order 



term in a series expansion of a solution to the Taub-NUT Einstein field equations (22 )-(25 ). From 
the multivaluedness, we see that u is ill-defined as a time variable, as we cannot access negative 
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values of u here, whereas P\{u) is multivalued as a function of u. Even more pathologically, P 2 
(which essentially plays the same role as g(u)) seems to originate only at u = 0, and then extend 
out to infinity in a straight line. In terms of the description in terms of time t , P 2 comes in from 
infinity, reaches the origin, and then traces back on itself. This already hints at an observation that 
will be made more precise below: our inability to access negative values of u presumably amounts 
to those values corresponding to some Euclidean, imaginary time coordinate obtained by a Wick 
rotation from the physical time t. Of course, the behaviour described here depends on the exact 
choice of the lower limit of the integration in the definition of u, but making a different choice 
does not avert the problem; it simply adds a constant offset to u. 

2. Direct solution of Einstein equations for reparameterised Taub-NUT metric 



Setting the scalar field F = (and A = 0) in (22)-(24), we can, in complete analogy with 



Section VII A[ solve directly for the Taub-NUT model parameterised in terms of Pi and P 2 . Note 



that we have already shown that for this coordinatisation there exists a simple scaling solution 



(26), contrary to the original Ryan-Shepley form (14) 
The vacuum equations 



AK[ + 8£ 2 - 2Kl + AK X K 2 + JL U^ - a\ = 0, (39) 



1*2 V '2 
2 



■> j-ua 

p}p} V p 2 



AK 2 + 6Ki + AK X K 2 - -„ 3-^-4=0, (40) 



1*2 

8 (P\ 



j2 



-2K 2 2 -A Ki K 2 + -^^-AJ=0, (41) 

can be immediately integrated by again taking an appropriate combination of a dynamical equation 
and the Friedmann equation to give 

Hu)= ( i±^i y, (42) 

where /i is a constant of integration. Note that P 2 (u) never goes to zero. In fact this form for 
P 2 (u) is nearly identical to the square-root of the standard form for g(u) in the Taub-NUT solution 



given in ( 14). The slight difference occurs because our time coordinate u differs to the one used 



by Ryan & Shepley, in general, by a relative shift uq. In (14), uq was chosen to be zero such 



that the function was symmetric around u = 0. Here, we will instead enforce our usual boundary 

21 




FIG. 8: Analytic Taub-NUT solution in terms of P x and P 2 for ji = /3 = 1. 



conditions on P\, namely that it passes through zero when u = 0, with a certain slope. Substituting 
the expression for P 2 into the Friedmann constraint allows us to integrate to find P 2 , with some 



additional integration constant /3. Now choosing uq 



J8(U 



such that Pi vanishes at u = 0, we find 



m») 



P 2 (u) 



p 2 li 2 u 2 -8pliu + p 2 + 16 J' 



p 2 li 2 u 2 -8piiu + p 2 + 16 
W 



(43) 



(44) 



Note the similarity with the derivation in Section VII A We could now eliminate the integration 
constants in favour of a§ and bo, or rather their equivalent series initial conditions. The functions 
P\ (u) and Piiu) are plotted for /1 = /3 = 1 in Fig. [8} 



3. Comparison of the vacuum solution in different set-ups 



It is clear from comparing ( 13 ) with ( 17 ) that g and Pf are in correspondence 



-P 2 {u)~g{u) 



and that the degenerate radii are straightforwardly related as 



1 



P 2 (u) ~ e 



C(«) 



(45) 



(46) 
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where the ~ is taken to denote equivalence up to the above shift uq. 



One sees that applying the above identifications (46) and (45 ) to the Taub-NUT scaling family 



(26 ), one obtains a family of solutions even in the conventional setup: 



g(u) = j=g(Pu), CW = C(/3")-^ln/3. (47) 



This can also be seen from changing the time parameter in the metric itself as before. However, 
this transform looks less natural than the one for DLH, adding to the suspicion that Taub-NUT 
time is not a physically sensible coordinate. 

One may continue to match the vacuum solution in the three different set-ups analytically. The 
above matching of the radii in the two versions of Taub-NUT together with the relative time shift 
uq allow one to identify A and B in terms of /i and /3. This now completes the identification of 
the two Taub-NUT versions. These two equivalent models can now in turn be related to the initial 
conditions qq and bo in the (original and periodic) DLH model using the time reparameterisation 



(15). 



In fact, all three models can be shown to coincide, at least for the portion where the radius 
in Taub-NUT corresponds to some physical distance (g(u) ~ Pf (w) > 0), i.e. where taking the 
square root gives something real. Fig. |9| displays yjg ~ P\ from their analytical forms ( 14) and 



(43 ) in the range where g(u) ~ Pf > ( y/g dashed, P\ dotted). They are found to be in agreement 



with our matching of these two solutions analytically above. Furthermore, they also coincide with 



a plot of R\{i) = P\(u) as a function of u(t) as defined by (15), evaluated for the periodic DLH 
vacuum solution considered above. However, whereas y/g and P\ only live in the upper half plane, 
the mapped DLH solution winds around a mirror-symmetric closed parametric curve in P\-u{t)- 
space as it cycles through its periodic oscillations, with one half-period of the elliptic solution 
corresponding to the positive portion of g. The other radii e% P2 and R2 match likewise. 

This now raises the question: to what do the parts of the analytic Taub-NUT solutions with 



g < correspond when mapped over to the DLH setup using the inverse mapping (21 )? Clearly, 



in this regime, the physical radius P\ = R[ must become imaginary. The prescription for finding 



t(u) then indicates that time t must likewise become imaginary, as follows. Using the relation (20) 



together with the known analytic solutions ( 14 ) and (43 ) we find 



A_^ = J> / j3W-8j3/x« + P 2 + 16 / J g(4gig+1) 

du P x y/g V J3w(-4J3jUm + j3 2 + 16) \ Au+l-4B 2 u 2 ' { ) 
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FIG. 9: Explicit matching between conventional Taub-NUT (y/g(u), dashed black), our reparameterised 
version of Taub-NUT (P\(u), dotted black) and the periodic vacuum DLH model mapped over using the 



diffeomorphism {R\(t) versus u{t) as given by (15), continuous gray). They are found to lie on top of 



each other, with the Taub-NUT models covering the upper half plane, corresponding to the region where 
g(u) ~ P\{u) > 0, whereas the DLH vacuum model winds around a closed curve in u space as it cycles 
through its periodic behaviour. The parameters are matched such that the periodic DLH model has the same 
initial conditions as considered earlier. The other radii e^ , P2 and R2 can be shown to match likewise. 

This is well-defined when g > and just recovers the positive half-period of the DLH vacuum 



model, as is obtained numerically (see Fig. 10). However, when g < we can perform a Wick 
rotation to imaginary time x = iu. In the limit of large u where the Taub-NUT solution approaches 
a negative constant, we find that t and u are linearly related, but with a relative factor of i, such 
that t ~ T ~ iu. So the NUT regions of Taub-NUT correspond to both imaginary space and time 
coordinates in the DLH model. Fig. [TT] displays the u < NUT region mapped over to DLH by 



integrating (48 ) numerically. Here, we have chosen to plot the result in the first quadrant, as there 
is a freedom of choosing factors of ±i on both axes. At late times, the linear relationship between 
t and u means that this mapped NUT region also settles down to a constant like in the original 
Taub-NUT model. So, in particular, mapping the Taub-NUT model does not recover the periodic 
DLH solution, and instead corresponds to imaginary space and time. Performing the integration in 



48]) analytically yields a solution in terms of elliptic integrals of precisely the same form as (33), 
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FIG. 10: The Taub region in Taub-NUT (corresponding to g > 0) mapped over to DLH space (black) falls 
on top of the the first half -period of the periodic vacuum DLH solution (grey). 



as one would expect. 

One marked difference between the DLH vacuum solution and this new version of the Ryan 
and Shepley form of the Taub-NUT solution is now that the latter is not periodic, which is the 
opposite of what one would commonly expect: normally Lorentzian models are not periodic, 
but upon euclideanising become periodic (cf. the periodicity in imaginary time of Green's and 
partition functions, in particular in relation to black hole thermodynamics). However, a similar 
phenomenon has recently also been observed in Bianchi V 11331 . 

Note that when one allows both the positive and the negative branches of Jg, one recovers 
precisely the trajectory in w-Pi-space that the mapped DLH model traces out (Fig. [9]). However, 
in the conventional setup Taub-NUT only selects one branch, and then turns to imaginary time 



and space on either side (the NUT regions, see Fig. 11). This suggests that even the Taub-NUT 
model could cycle indefinitely (like the vacuum DLH model) if it was allowed to change from 
one branch to the other. In fact, this matches smoothly: we have chosen Pf to have a zero at the 



origin, so it is linear there from (43) . Thus Pi goes like -Ju at the origin, as also observed in 
Sections VII C 1 (Fig. [7]) and VII D| Hence, there is actually a rather natural, smooth transition 
between the two branches, exactly like in Fig. [7] This also suggests that the smooth, single- valued 
parabolic behaviour of DLH is more natural and more fundamental, as opposed to the Taub-NUT 
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FIG. 1 1 : First NUT region in Taub-NUT (corresponding to u < 0) mapped over to DLH space 

case, where the parabola essentially gets turned sideways such that the description in terms of 
square roots results in multivaluedness by having the two different branches. 

D. Taub-NUT with a scalar field 



Now that we have clarified the relationship between the DLH and reparameterised Taub-NUT 



set-ups in the vacuum case, we can use the transformation (15) to map our cosmological DLH 



solution with a scalar field, presented in in [1J, directly over to the reparameterised Taub-NUT 



model. Alternatively we could find an analogous series expansion to (|22|)-(|24|) directly. 
In either case, to lowest order, the series solution has the form 

P\(u) = y/u, Piiu) = const, F(u) = const, 



(49) 



and thus analogous boundary conditions to those used in the DLH case may be used, from which 
a numerical integration can be performed straightforwardly. 

Gradually introducing a light or diffuse scalar field has a similar effect to the analogous scenario 



in DLH (Section VII B), where it slightly perturbs the vacuum case. Fig. 12 shows a plot of an 
interesting choice of initial parameters with a sizable scalar field (arj = 1.0, bo = 17, K = 1, A = 0, 
fo = 6, m = |), that exhibits inflation and isotropisation thereby looking rather like the original 
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FIG. 12: Taub-NUT inflation: inflation and isotropisation in analogy to the DLH scenario for the model 
with parameters a® = 1, bo = 17, K = 1, A = 0, /o = 6, m = I. 
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FIG. 13: Taub-NUT inflation: inflation, but not complete isotropisation for the model with parameters 
ao = 1, bo = 2, K = 1, A = 0, fo = 6, m = 1. 
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DLH solution. Fig. 13 shows a plot of another interesting model with the set of initial parameters 
ciq = 1.0, bo = 2.0, K = 1, A = 0, fo = 6, m = 1), that exhibits inflation, but does not completely 
isotropise. However, the different Hubble factors tend to a common value so in the flat-space 
late-time limit the difference in scale factors would be unobservable as argued in [Q]|. 



The slope of P( in Fig. 12 at the pancake is unity, so that slope Pi here is 1/2, in agreement 



with the DLH series solution mapped to Taub-NUT at this point (see equations (49) and (38)) 



The solution in Fig. 13 does not appear to be of the pancaking type, so we do not consider its 
early slope here. However, the late-time slopes in P\ are 0.4 in both models (i.e. 0.8 in P% ). Note, 

2 

that this does not agree with the behaviour a(t) oc fi that one would expect once the scalar field 
has decayed and behaves like non-relativistic dust in an Einstein-de-Sitter phase (at which point, 
presumably, reheating would occur). Note, however, that the DLH model in Fig. [Tjdoes have the 
required late time slope of |, so that the DLH coordinate time matches onto physical time at late 
times. So we conclude that the Taub-NUT time coordinate is not a physical time coordinate at late 



times, unlike the DLH time. This is consistent with the coordinate transformation ( 15 ). 



We can compare the slopes in the DLH model and the corresponding Taub-NUT model for 
late-time power-law behaviour R(t) = at? as follows. The corresponding slope in Taub-NUT is 



(50) 



dP 


dRdt 


dR2 


2/3 


du 


dt du 


' dt R~ 


t 



When we can assume that the integral in ( 15 ) is in fact dominated by the contribution from the 
power-law behaviour we can integrate to get 



a 



2(0 + 1) 



tP+i 



Substituting (51 ) into (50) then yields 






+ r 



2Q6 + 1) \J~r 

11 

a 



"|8+1 



(51) 



Integrating we obtain 



du \ 



P{u) = a^ (2(JS + 1)) P+ 1 WP+ 1 +c. 



(52) 



(53) 



For an Einstein-de-Sitter phase with /3 = I this therefore gives P{u) ~ u 2 ' 5 , which is indeed what 



is observed in Figs. 12 and 13 



Fig. 14 shows u(t) according to ( 15 ) in both the original DLH model with a scalar field and 
the new periodic vacuum solution. Note that the early time slope in u(t) is 2, in agreement with 
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FIG. 14: u(t) for the original DLH model with a scalar field (left) and the periodic vacuum DLH model 
(right). 



formula (51 ) in the DLH model. The slope at late times in the DLH model is 1.67 ~ 1 + |, also 



in agreement with (51 ). Note that the integral ( 15 ) is indeed dominated by the contribution after 
inflation. As already mentioned in Section |v| in the periodic model u{t) can be seen to go back 
on itself when R\ changes sign. This means that u does not measure the progression of time 
in the same way as t, but instead retraces a parameter range that has been traversed previously. 
Considering that DLH time t matches onto physical time, this casts further strong doubts on Taub- 
NUT 'time' u being a sensible measure of time. 

We conclude that the desirable features of inflation and isotropisation can survive the mapping 
to Taub-NUT. However, at early times - near the pancaking events - the DLH coordinatisation 
seems much more natural (c.f. Figs [7]). At late times, the physical time coordinate can be inferred 
from the expansion history of the universe, and Taub-NUT again fails to match onto something 
physical. Furthermore the scaling solution is most natural in DLH, and very contrived in Taub- 
NUT coordinates. This also ties in with the fact that the closed Bianchi model is an immediate 
generalisation of the physical closed FRW model, which also admits such a scaling family. The 
multivaluedness of the Taub-NUT time coordinate when mapped from a physically meaningful 
model is a problem, as is allowing it to go imaginary in order to reconcile the different solutions. 
Thus Taub-NUT time u fails to be a physically sensible measure of time on several accounts, and 
we therefore contend that the DLH coordinates are the more physical coordinate system. 
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VIII. CONCLUSIONS 

We have compared and contrasted two biaxial Bianchi IX spacetimes - Taub-NUT and the DLH 
model. They exhibit great similarities such as dimensional reduction at a non-singular pancaking 
event, inspiralling geodesies and the same eigenvalue structure of the curvature and Weyl tensors. 
However, there are profound differences in the global structure and physical interpretation. We 
have shown that the two metrics can be mapped into each other, but that the coordinate transforma- 
tion is multivalued. This property is responsible for introducing artifacts into the coordinatisation 
and thus accounts for the differences in global behaviour. We believe that our parameterisation 
in terms of physical scale factors in complete analogy with, for instance, FRW-models is more 
natural, as opposed to working with the squares thereof and allowing those to become negative. 
In the light of this, we have removed the scalar field from the original DLH model and found an 
analytic vacuum solution in terms of elliptic integrals. This solution is periodic and also recovered 
by numerical integration, with which there is very good agreement. There is also an obvious link 
between this periodic solution and the pancaking series solution of the DLH model, in that the 
boundary conditions at the pancake are the same. However, after the inclusion of the scalar field 
the DLH model isotropised, inflated and yielded a viable late-time cosmology with a perturba- 
tion spectrum consistent with observations as shown in JTJ]. Thus in either case, working with the 
scale factors yields sensible cosmological models and, in particular, the spatial sections are closed 
Bianchi IX throughout, in contrast with the topology changing transition in Taub-NUT. Mapping 
the periodic DLH vacuum solution using the coordinate transformation matches onto the analytic 
Taub-NUT solution for the Taub portion where g > 0. However, in the NUT regions, where g < 0, 
we find that Taub-NUT time corresponds to imaginary time in 'DLH space' . This firstly adds to 
the doubts that it is a sensible time coordinate; secondly, it resolves the alleged topology transi- 
tion, as DLH now has closed spatial sections throughout as well as a physical time coordinate. 
Thirdly, and rather surprisingly, this suggests that the Lorentzian Taub-NUT is actually Euclidean 
in these regions. It is then surprising that the (Lorentzian) DLH model is periodic, whereas the 
'Euclidean Taub-NUT' is not, which is contrary to what usually happens. We also note that DLH 
time essentially acts as conformal time for Taub-NUT. 

The DLH coordinatisation thus has a number of advantages over traditional approaches. Firstly, 
the behaviour of the solution across pancaking events seems much more natural and straightfor- 
ward in both the DLH vacuum and scalar field case. The well-behaved odd-parity series solution 
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that we had found previously acquired bad behaviour when mapped over to Taub-NUT due to the 
fact that the transformation is singular at pancaking events. Secondly, the periodicity of the DLH 
vacuum solution appears more fundamental than the Nut-Taub-Nut structure and in fact it is clear 
how this periodicity gets broken by an integral involved in the mapping. Thirdly, it is clear that 
the scale factor has physical significance and its square should therefore not take negative values. 
The resulting alternative Taub-NUT parameterisation still has the advantage of admitting a simple 
scaling relation over the conventional metric, which only admits a rather awkward looking scaling 
family. Moreover, at late times when we can infer cosmological time from physical measurements, 
the DLH model behaves as one would expect physically, whereas Taub-NUT time again fails to 
produce physical results. We conclude that our coordinatisation provides at least an interesting 
alternative point of view that could shed some light on some of the pathologies that Taub-NUT is 
believed to have, and might actually be the more natural coordinatisation. 
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Appendix A: Mapping of curvature invariants 

Given the similarities in the metric Ansatz, i.e. biaxial Bianchi IX on the hypersurfaces, it was 
obvious that similarities between the Taub-NUT and DLH geometries should be explored. In fact, 
there does exist a coordinate transformation linking the two spacetimes, which we have presented 
above. In general, however, when one wishes to examine the relationship of two spacetimes it is 
often the easier strategy to show that the two spacetimes are actually distinct. This is made difficult 
by the general coordinate invariance of General Relativity, so one must find coordinate invariant 
ways to distinguish between spacetimes, such as curvature invariants. A convenient classification 
is the Petrov classification, which examines the algebraic properties of the Weyl tensor. This is 
often easier than other approaches due to the additional self-duality structure and tracelessness of 
the Weyl tensor (see, for instance, [|34ll for a reference). If, however, it turns out that the Petrov 
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classification is not enough to distinguish between the two spacetimes (here both are of type D) 
one must find more curvature invariants, such as topological invariants, or principal curvatures. If 
after this kind of tests one has still failed to demonstrate that the two spacetimes one wishes to 
compare are distinct, then the curvature invariants can be used to find the diffeomorphism linking 
the two spacetimes, if indeed such a diffeomorphism exists. 

The Riemann and Weyl tensors R^vat and C^vctt are multilinear operators of fourth rank acting 
on tangent vectors. However, they can also be considered as linear operators acting on bivectors, 
and as such they have a characteristic polynomial, whose coefficients and roots (eigenvalues) are 
polynomial scalar invariants. That is, consider the eigenbivector equations 

"fiva'c^ = ^"J/jvi Levari = i*-\lvi (Al) 

for some bivectors 5^ v and T^ v and eigenvalues A and 7 respectively. The degeneracy structure of 
the eigenvalues provides a coordinate-invariant way of distinguishing spacetimes. Consideration 
of the eigenbivector equation of the Weyl tensor leads to the Petrov classification, which is in 
general easier to examine due to the additional self-duality structure of the Weyl tensor. The 
eigenvalues of the Riemann tensor are also known as the 'principal curvatures'. 

1. Petrov type 

There are six bivectors in 3 + 1 dimensions, so that, in general, a fourth rank tensor acting on 
them has six real eigenvalues. Alternatively, the self-duality of the Weyl tensor yields a natural 
complex structure, so that the six real eigenvalues are equivalent to three complex eigenvalues. It 
turns out that both DLH and Taub-NUT have one degenerate pair of complex eigenvalues, which 
then also fixes the remaining one due to the tracelessness of the Weyl tensor. This corresponds to 
Petrov type D, which is what Taub-NUT was previously known to be. In terms of real eigenvalues, 
there are two degenerate pairs and two singlets. The complex structure also results in the pairs 
being complex conjugates of each other, and likewise for the singlets. The real eigenvalues are 
found to be of the form 

yf = yD = Cl + v /Q, y° = yj> = C i -y/C 2 , $ = C 3 + VQ, ^ = C 3 -v^4, (A2) 

for the DLH model, with some expressions C, in terms of the 7? ; s and their derivatives. There is 
one constraint among them since the Weyl tensor is traceless, but we have suppressed the exact 
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results in order to avoid unnecessary clutter. Analogous results hold for the Taub-NUT model for 
some constants Z),, from whence the eigenvalue structure (Petrov D) and tracelessness can be seen, 
and the two sets of results map into each other as suggested. 

Alternatively, one can consider the Weyl tensor as a function W acting on bivectors B, such 
that the eigenvalue equation is W(B) = yB. The self-duality translates into W(IB) = IW(B) for the 
complex structure denoted by I. In terms of complex eigenvalues, the eigenvalue equation becomes 
W(Bi) = ((Xi + PiI)B(. Solving this in the DLH and Taub-NUT setups yields the eigenvalues 

a? = -2a? = -2af , al = -2a\ = -2aJ, j8 3 D = -2jB? = -2j8f, /3 3 r = -2j8f = -2$, 

(A3) 
which also enjoy the appropriate degeneracy structure, mapping properties and tracelessness. 

2. Principal curvatures 

The Petrov classification might not be sufficient to distinguish between two spacetimes, in 
which case one can consider the eigenvalues of the full Riemann tensor, which are also called 
'principal curvatures'. These are more complicated, as the self-duality (and therefore natural com- 
plex structure) and tracelessness are lost. However, it again turns out that both models have two 
degenerate pairs of real eigenvalues, and two singlets, mirroring the structure of the Weyl tensor. 
Computation of the eigenvalues gives the following form 



X l = A 3 = E\ + v^2, A 2 = A 4 —E\ — v£2, A 5 = £3 + \/Ea n A 6 = £3 — VE4, (A4) 

for the DLH model, again with analogous results for the Taub-NUT case. It can be seen that these 
share the same degeneracy structure of the eigenvalues, which is similar to the one of the Weyl 
tensor. The two models can also be explicitly mapped into each other, as suggested. 

Appendix B: Mapping of geodesies 

In order to complete the mapping between DLH and Taub-NUT, we now turn to showing the 
equivalence of the geodesic equations in both models. The geodesic equations are most easily 
obtained using the Lagrangian formalism, in which 

L = 4g^x v (Bl) 
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is varied with respect to the coordinates [jc /x ] = (t,x,y,z)', here a dot denotes a derivative with 
respect to some affine parameter X and the factor of 4 is included for later convenience. An 
explicit realisation of the Maurer-Cartan forms is, for instance, 

ft) 1 = dx + sinydz, 

CO 2 = cosxdy — sinxcosydz, 

ft) = sinxdy + cosxcosydz- (B2) 

Inserting the DLH metric ([6]) into L yields 

L D = 4i 2 - R 2 x 2 - 2R 2 smyxz-RJy 2 -(R 2 sin 2 y + R 2 , cos 2 y)z 2 . (B3) 

Since L is independent of the x- and z-coordinates, the corresponding Euler-Lagrange equations 
yield two conserved quantities K x and K z according to the relations 

-2R 2 (x + sinyz) = K x , (B4) 

-2[R 2 l smy(x + smyz)+R 2 ,cos 2 yz] = K z, (B5) 

which may be solved for x and z to yield 

K Z R\ sin y - K x (R 2 sin 2 y + R\ cos 2 y) 
2i? 2 i? 2 cos 2 y ' 

J^siny-JT Z 



(B6) 



L 



9P2 2 • ( B7 ) 

ZK2COS z y 

Substituting these expressions back into the Lagrangian, we get the modified form 

D _ 1 16f 2 i? 2 ^cos 2 y-4y 2 i? 2 ^cos 2 y- (K x 2 + K 2 ) R\ + 2R\smyK x K z + K 2 (R\ -RQ cos 2 y 

~4 i? 2 i? 2 cos 2 y 

(B8) 

The Euler-Lagrange equation for y then reads 

AR\ cos 3 y \R 2 y + 2 -^fy J = {K 2 + K 2 ) siny - 2K X K Z + K X K Z cos 2 y, (B9) 

as stated and further analysed in JTJ. 

For Taub-NUT we keep the same variables for now, except t is replaced by u. (In general, the 1- 
forms would be of the same form, but the coordinates and the affine parameter could be different.) 
We will see later that this choice is indeed consistent. 

L = A gpLV x^x v (BIO) 
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is varied with respect to the coordinates [x* 1 ] = (u,x,y,z). Inserting the Taub-NUT metric in the 



form ( 17 ) into L yields 



U = $u(x + sin yz) - P 2 x 2 - 2P{ sin viz - P^f - (if sin 2 v + P 2 cos 2 y) z . 



(Bll) 



Since this Langrangian is likewise independent of the x- and z-coordinates, we still have the 
two conserved quantities K x and K z . However, they pick up an extra term in u 

Su-2P 2 (x + sinyz) = K x , (B12) 

8wsiny--2 [Pf siny (x + sinyz) + P% cos 2 yz] = K z = K x siny-2P 2 cos 2 yz, (B13) 

which may be solved for x and z to yield 

Au K z Pl sin y - K x (Pf sin 2 y + P% cos 2 y) 



K x s'my — K z 



2P 2 P 2 



'f/^2 cos y 



(B14) 
(B15) 



2i , 2COS 2 J 

Substituting these expressions back into the Lagrangian, we get the modified form 
T 1 64ii 2 J P 2 cos 2 }; - 43) 2 J P 1 2 / 5 2 4 cos 2 ); - (K 2 +K 2 ) P 2 + 2P\ sin yK x K z + K x 2 (P 2 -P%) cos 2 y 
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P[P% cos 2 y 



(B16) 

where the terms that are linear in u have cancelled, and the use of the conserved quantities has in 
fact introduced a term that is quadratic in u. 



Comparing with (B8), we see that a straightforward identification of coordinates, scale factors 
and affine parameters is indeed possible, and the differing terms 64 ji 2 /^ cos 2 y and 16 t 2 R 2 R 2 cos 2 y 



are consistent with the time rescaling (20). Therefore the equivalence of the geodesies is already 



exhibited at the level of the Lagrangians. The Euler-Lagrange equation for y in Taub-NUT is in- 



deed identical to ( |B9[ ), as the Lagrangians only differ on the time-time piece, which is independent 
of y and y. 



Appendix C: Derivation of the elliptic integral solution for the vacuum DLH model 

In order to recover the analytic form of the numerical solution in Fig. [3] we note that an appro- 
priate combination of (29) and (31 ) allows us to integrate a Bernoulli-type equation in Ri(t)~ 2 to 



find R{ in terms of Rj and an integration constant ji as 

R2(t)R 2 (t) 



*i(0 



~y/nR 2 (t) 2 -l 

35 



(CI) 



As will become clear momentarily, we will choose the minus sign, in order to recover the analytic 
form of the numerical solution above. Substituting this expression for R\ into equation ( f29| ) allows 
us in turn to solve for R 2 . This gives another integral solution with positive and negative branches 
and two integration constants. Differentiating and choosing the branch that is decreasing at t = 0, 
we lose one integration constant and can fix the other in terms of /! and bo such that R 2 (0) vanishes 
in line with the pancaking series solution and the numerical solution above. From ( |C1[ ), R\ will 
then be increasing and vanish at t = as asserted earlier. The slope of R\, which is by definition 
ao, can then be used to fix \x and hence the whole solution in terms of initial conditions 

an + 4 
M = J W- (C2) 



ci 



2 

0^0 



b 2 ' 



With the constants in the problem explicitly fixed, we now return to the integration of R2 and R\. 
From the denominator in dCT} this will involve finding integrals of the type 



'a 2 R 2 (t) 2 + 4R 2 (t) 2 -a 2 b 2 = x. (C3) 

In fact it can be easily seen from the numerical solution that bo simply scales both axes, so we set it 



to unity in the following. Inverting (C3 ) and substituting in for R 2 in terms of x gives a dynamical 



equation for x that can be integrated to find time t in terms of x (and thus R 2 ) as 



' - X ^ a K+xf^E(y;i) + x f^(ao + 2i)F(yU)-^(al-4)n(y-,^-,i). (CI) 



2an a.Q + 4 \ Si ' V 32 ' y 32an \ ' ^0 

where the elliptic integrals have been reduced in terms of Legendre's three normal forms 



F( * i) = f* ^-M-eA <C5) 

n(z;v;£)= / F = ds. 



Implicit in the Legendre forms is a choice of branches. In ( |C4| ) these are chosen such that 
Maple's default choice indeed recovers a real solution as required. The explicit numerical eval- 
uation of this analytic solution in Maple is confirmed to be real (to within numerical precision) 



and Fig. 15 shows very good agreement between the analytical solution and the results from the 



numerical integration for four values of ao: ao = 0.5, 1.2, 2, 5. Time t has been calculated as 
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#2 



ao = 5 



a = 2 




ao= 1.2 



«o = 0.5 



FIG. 15: Comparison of analytic and numerical solutions (which coincide) for 4 values of ao'. ao = 
0.5, 1.2, 2, 5. For each of these values the analytic expression is indeed real and coincides with the nu- 



merical result to great accuracy. Note that due to the critical factors of y/ao — 2 in (C4i there can be relative 
signs depending on which branch one chooses in the function definitions as one crosses «o = 2. Here, the 
last two terms have flipped signs for ao < 2 when explicitly evaluated (using Maple) relative to the analyti- 



cal ones in (C4i. Time t is calculated as a function of /?2 in the region between the initial maximum and the 
first minimum of R2 due to the multivaluedness and then plotted along the horizontal for comparison with 
the numerical results. 

a function of R2 between the first maximum and minimum of R2 due to its multivaluedness, and 
then plotted along the horizontal axis for comparison with the numerical results. 

If one repeats the above analysis for bo 7^ 1 one finds that bo always occurs with the same 
power as x as might be expected from the definition. Changes in bo can therefore be absorbed into 
x. Moreover, the arguments in the elliptic functions turn out to have overall power whereas the 
coefficients have weight one. This explains how the time t gets simply rescaled. Of course, R2 
also gets rescaled, given that bo is its boundary condition at t = (which is effectively the same as 
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absorbing it into x). 
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